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We consider the possibility of detecting spin precession in a magnetic field by nonequilibrium 
transport processes. We find tfiat time reversal symmetry imposes strong constraints on the problem. 
Suppose the tunneling occurs directly between systems at two different chemical potentials, rather 
than sequentially via a third system at an intermediate chemical potential. Then, unless the magnetic 
fields are extremely strong or spin polarized electrons are used, the periodic signal in the current 
results from beating together two different precession frequencies, so that observing a signal near 
the Larmor frequency in this case requires having some cluster with a g factor close to zero. 



PACS Numbers: 03.65.Xp, 73. 23.-b, 73.63. Kv 

Transport of spin-polarized current through an STM 
[1,2] opens the possibihty of investigating magnetic sys- 
tems at a much smaller scale than is possible using elec- 
tron spin resonance [3] . The ability to manipulate single 
spins via STM is a key part of the development of spin- 
tronics [4] and quantum information processing. 

Interestingly, even using spin-unpolarized electrons, 
a periodic modulation of tunneling current through an 
STM tip was observed experimentally [5]. The original 
proposed explanation depended on periodic modulation 
of charge at the Larmor frequency by a precessing spin 
center. Another very interesting possibility depends on 
spin-orbit scattering [6-8]. However, in this paper, we 
will show that, with two exceptions, time reversal sym- 
metry imposes strong constraints on the problem, and for 
both of the mechanisms considered leads to modulation 
at a frequency resulting from beating together two dif- 
ferent precession frequencies. Then, for a single impurity 
spin, the modulation in the current appears at twice the 
Larmor frequency, and not at the Larmor frequency it- 
self. To obtain a signal at the Larmor frequency requires 
two spins, one of which has a vanishing g-factor [9]. 

The first possible exception is to have a sufficiently 
large magnetic field to introduce Aharonov-Bohm phases 
for the electrons; in this case time-reversal symmetry of 
the tunneling Hamiltonian is broken by the external field. 
This case is unlikely to be experimentally relevant, due to 
the strength of fields required. The more interesting ex- 
ception involves sequential tunneling through a quantum 
dot [6,10], and will be discussed more below. Sequential 
tunneling refers to a system in which there are two leads 
and a quantum dot, with the chemical potential one lead 
far above that of on the dot, and the chemical poten- 
tial on the other lead far below that of the dot. Thus, 
this system has three different chemical potentials; the 
results in this paper will apply to systems in which there 
are only two different chemical potentials. 

As a model, we first consider a tunnel junction with a 
single nearby spin. The electrons can couple to the spin 



via exchange coupling, and there may also be a spin-orbit 
coupling present. This provides a means of measuring the 
spin in terms of its effect on the tunnel current. We first 
consider the appropriate Hamiltonian for the system in 
the absence of any external magnetic field, so that the 
system is time reversal invariant, and then consider the 
effect of any applied magnetic field. We find that the ef- 
fect of time reversal symmetry strongly limits the allowed 
modulation of the current by the spin. Our physical setup 
with an impurity spin and a single tunnel junction will 
be similar to that considered in [13]; similar results with 
some exceptions hold for tunneling through a quantum 
dot as considered in [6], as discussed below. 

Time Reversal Symmetry — Consider the most general 
possible Hamiltonian, written as H = "Hl + "Hr + T^t, 
where TLl.r are the Hamiltonians for the left and right 
leads of the system and Tir includes the tunneling terms, 
as well as terms coupling the electrons to the spin. The 
two leads may be taken to be Luttinger liquids. There is 
a voltage difference V between the leads, giving a chem- 
ical potential difference /i = qV , where q is the electron 
charge. We assume that the effects of electron interaction 
are negligible in the tunneling region, as the tunneling 
process is weak and only one electron tunnels at a time. 
Thus, we restrict to terms which are at most second order 
in the Fermion operators in TLt- 

Then, given the requirement of time reversal symme- 
try, the most general possible Tir, is 



TLt = + tiJ^Vj + t2-h{SjWj) 



(1) 



where a sum over repeated indices is assumed. Here 
Sj (for j — X, y, z) is the spin operator on the impu- 
rity spin (which is taken to have total spin 5, where S 
may be larger than 1/2) and S^^"^ = a^'^Sj, where aj 
are the Pauli spin matrices. The letters /3, 7 identify 
the particular lead: — L,R, while ip'^ are electron 
destruction operators with spin fi =t,i- We have de- 
fined the following operators: first, a tunneling opera- 
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tor To = {''Pi)'^ '4>'ii + h.c. Next, an electron spin current 
Jq = i(Tj'^[{ijjjj)^ip'^ — h.c.]. Finally, an operator related 
to the current operator, Jq = i[{tp'^)''tp'^ — h.c.]. The term 
in to is a tunneling term. The term in ti, which combines 
spin-flip with tunneling, can arise via a mechanism sim- 
ilar to that generating spin scattering by nonmagnetic 
impurities in semiconductors [12], as noted in [6]. The 
fixed vector Vj is determined by the spin-orbit scatter- 
ing. The 4 couplings J^^ are exchange couplings, while 
the fixed vector wj is set by both spin-orbit and exchange 
interactions. 

The terms in (1) are time reversal symmetric [11]. This 
operation complex conjugates all terms in the Hamilto- 
nian, and sends ip'^ —>■ and —>■ —tp\ and similarly 
for ip^. Further, it changes the sign of the spin operator. 
Si —Si. Then, Jq changes sign under time reversal as 
does S, while Jq is time reversal invariant. 

The coupling of the tunneling current to the spin de- 
pends on a combination of exchange interaction between 
the electrons and spin with spin-orbit scattering. Since 
both of these interactions are in the term with t2, we 
focus on the following Hamiltonian: 

(2) 



We = + 'Hr + toTo^ + t2JoSz 



Here, we have chosen w to lie along the z-direction. How- 
ever, many of the results we obtain will be valid for much 
more general Hamiltonians, as will be discussed at the 
end. 

External Magnetic Fields — To Eq. (2), we add an ex- 
ternal magnetic field. This field breaks the time reversal 
symmetry and allows for terms which are odd under time 
reversal, so that Ji = He + Ti-o- Two particularly impor- 
tant allowed terms are 



Wo = BiSi + tsJo- 



(3) 



The first term is the coupling of the field to the spin, while 
the second results from a possible Aharonov-Bohm phase 
on the electrons traversing the tunnel junction: the phase 
of the term coupling to Sz has been changed by the mag- 
netic field, so that it is no longer imaginary. We expect 
that terms in t^ can be disregarded for experimentally 
relevant fields. 

The possible introduction of terms like is very im- 
portant. To see this importance, return to the Hamil- 
tonian Tig. In this Hamiltonian, the only spin operator 
that appears is the operator Sz. Thus, if we pick a basis 
for 5* that diagonalizes Sz we can consider S to be con- 
stant in time, and thus Sz in the Hamiltonian (2) can be 
considered as a classical parameter. Then, this Hamil- 
tonian is equal to Hl + Hr + [^(V'l)^^^ + h.c], where 
z = to + it2Sz. Since the phase of z is unimportant, phys- 
ical parameters like the current can depend only on the 
modulus: j^j = \to + it2Sz\ =tl + t^S"^. 

Thus, the current depends on 5^, but cannot depend 
on the sign of Sz. In fact, for the most general Hamil- 
tonian of Eq. (1), this remains true, as will be discussed 



more below. On the contrary, if we consider the com- 
bined Hamiltonian, He + t^Jo, the modulus becomes 
t'o + t%+2t2t3Sz+t'2Sl, and thus does depend on the sign 
of Sz ■ This difference in the dependence of the current on 
Sz, when Sz is considered as a classical field (unchanging 
in time) determines, as we will see, how strong the modu- 
lation of the current is at the Larmor frequency, ujl , when 
Sz acquires time dynamics due to the introduction of an 
external magnetic field. The Hamiltonian He+Ho, with 

7^ 0, has been considered previously in the context of 
a two level system rather than a spin, in which case such 
terms do not violate time reversal symmetry [13], as well 
as in the context of a spin system for which such terms 
do violate time reversal symmetry [7,8]. These authors 
found a modulation in the tunneling current at the Lar- 
mor frequency. In the case of tunneling through a dot, 
as considered in [6], similar time reversal symmetry vio- 
lating terms were considered in the matrix O^ss' of that 
paper, again finding modulation at the Larmor frequency. 

Current Modulation at Weak Field — We now consider 
the case instead with a sufficiently weak field, so that t^ 
can be disregarded: H = He + BiSi. We will find only 
a very weak modulation of the current at the Larmor 
frequency, and a much more substantial modulation at 
twice the Larmor frequency. This has important experi- 
mental implications for the detection of a precessing spin 
with a tunneling current: for small tz, a much stronger 
signal should be found at twice the Larmor frequency in 
this setup with a single tunnel junction. 

The effective dynamics of the density matrix for the 
spin can be found using previous results [15], based 
on a Keldysh technique [16]. Throughout, we work 
perturbatively in the tunneling, Ht- Define /(/i) = 
27rpipi^/i"+^/r(a + 2). The parameter a depends on 
the Luttinger parameters of the leads and the sample ge- 
ometry; it is equal to zero for Fermi liquid leads [14]. The 
product of the density of states in the left and right leads 
at energy fi is plPrI^" ■ Then, the result for the effective 
dynamics is 



t2 d/ 
2 d/i 



,{a,p}] 



t 



s. 



(4) 



where A = [Hq, t2Sz], and Ho is the Hamiltonian for the 
spin. This includes the external field, BSx, as well as an 
additional effective field along the z direction due to the 
coupling to the electron current. We choose to absorb the 
additional field into a renormalization of the bare field, 
so that the net field is along the x-direction: Ho = BSx- 
Thus, A = -it2BSy. 

The density matrix evolution (4) describes the Hamil- 
tonian evolution of the spin under Ho, as well as dissipa- 
tion (the second term) and decoherence (the third term). 
For large V, Eq. (4) leads to an effective temperature 
Tefi — {a + l)~^qV/2. This gives a stationary distribu- 
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tion of the density matrix for the spin: 

P=^^(l-5^x/Teff). (5) 

Current — We now consider the average current across 
the system. The current operator J is defined by taking a 
derivative of the Hamiltonian with respect to an external 
gauge field, giving 

J = q{toJo-t2ToS,). (6) 

The average current is, for B « fi, [15] 

(J) = q{tl + tlSl)I{,i) - qBtl{[S,, (7) 

The first term is the current that results from consid- 
ering Sz to be a classical object, while the second term 
is a weak quantum correction that slightly reduces the 
current due to fluctuations in Sz. 

As long as the net magnetic field is along the x- 
direction, (S'^)=0, and {S^) = (1/3)5(5 + 1). If the 
field acquires a 2;-axis component, then, for S > 1/2, the 
expectation value {SD will increase, and thus the net cur- 
rent will increase. The expectation value {[Sz, Sx\Sz) = 
(B/Teff)(l/6)S'(S' + 1). Thus, for a net field along the 
a;-direction, the net current is 

{J)=q{tl + {l/2.)S{S + l))I{n) (8) 
-qti{B^/T,^){l/&)S{S+l)^^. 

Current Noise — We now wish to compute the current- 
current correlator, {J{lo)J{—lo)). At second order in 
we obtain the shot noise contribution. It is ob- 
tained by using Eq. (6) for J and computing the cor- 
relation function {J{uj)J{—uj)) using the Hamiltonian 
Ti-L + Ti-R, and no higher powers of Ht- The result is 
(J(a;)J(— w)) = 2q{J), for « M- 

At fourth order in Ht, many contributions are possi- 
ble. The full calculation is quite involved, however the 
most important contribution is a "modulation" of the 
current. We consider a modulation of the current to be a 
term in ( J(a;) J(— w)} which is of order In con- 

trast, for example, the shot noise is only of order 

Such modulation terms are obtained as follows [15]: 
we first consider Sz to be a classical field, and ob- 
tain the dependence of the current on Sz, giving J = 
q{tQ + t2Sz)I{n)- Then, we consider how Sz varies in 
time, using the dynamics (4), to compute the correlation 
functions of Sz- The result is 

{JiO)Jit)) =qHili„f{S^MSm- (9) 

For S > 1/2, this gives a modulation of the current at 
twice the Larmor frequency. However, there is no modu- 
lation at the Larmor frequency itself for any S. 



Diagrammatically, we show a term giving such modu- 
lation in Fig. 1(a). Lines with arrows represent electron 
propagators, while pairs of straight lines without arrows 
represent spin propagators (here, simply to give some 
diagrammatic way of writing interactions with spins, we 
have chosen to use the standard represention of a spin op- 
erator as a pair of fermion operators) . The curly lines the 
current operator, while the crosses represent the operator 
Ht; since these two operators have straight lines leaving 
them, this indicates that we take the term —qt2ToSz in 
the current operator, and the term t2JoSz in Ht- There 
are two pairs of lines connecting the bubbles, giving a 
modulation at twice the Larmor frequency. In Fig. 1(b), 
we show a similar diagram with only one pairs of lines 
connecting the bubbles; however, this diagram vanishes. 

The fact that the modulation is not present at fre- 
quency (jJl will be true to all orders. In contrast, other 
terms in {J{—co)J{u)), which are lower order in /x and 
thus do not represent a modulation for this experimen- 
tal setup, may show a peak at the Larmor frequency. 
These terms all depend on Sz having a time dependence 
on time scales of order /i^^ (the time over which a tun- 
neling event occurs), and thus all of these terms will be 
suppressed by powers of B/^ and are therefore expected 
to be qinte small, for this situation. In such terms, we 
would find diagrams with both current operators on the 
same electron bubble, rather than on different bubbles as 
was the case in Fig. 1. 

However, for other situations, these terms may be im- 
portant. Consider sequential tunneling from one lead to 
the other via a quantum dot [6] , where the chemical po- 
tential in one lead is far higher in energy than the dot, 
and the chemical potential in the other lead is far lower 
in energy than the dot. In this case, the current through 
the system is largely independent of the exact difference, 
fx, between the chemical potentials in the two leads. Our 
considerations in this section relied on computing the or- 
der of terms in /x. However, in the given experimental 
setup with the quantum dot, the current is of order 
and the time to tunnel from left to right k^aci is deter- 
mined not by fi, but by the hopping matrix elements be- 
tween the quantum dot and the leads. Then, it has been 
shown that [10] for tunneling through a quantum dot in 
this sequential tunneling regime, one can obtain a con- 
tribution to the current at the Larmor frequency which 
is comparable to the zero frequency component of the 
current. Thus, while in [6] terms were considered which 
violate time reversal symmetry, with slight modifications 
a signal at the Larmor frequency can be found using a 
Hamiltonian which respects time reversal symmetry. 

More General Hamiltonians- The reason for a lack of 
modulation is that if the impurity spin is held fixed, 
with a static Sz, then the current depends only on S^ 
and not on Sz itself. For the Hamiltonians considered 
above, this results from the fact that the current de- 
pends only on the modulus of the hopping coefficient 
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z above. However, this is true even in a more general 
case. Consider any Hamiltonian including potential, in- 
teraction, and position-dependent coupling of the elec- 
tron current to S^,: U = i^^x){d'^ /2m. + U{x))tp{x) + 
^^{x)i^{x)i^^{y)^{y)V{x-y)+i^l^^{x){F{x),d,}^{x)S,. 
The spin-orbit coupling introduces an effective vector po- 
tential. If this vector potential has a non- vanishing curl, 
then it leads to a non-vanishing effective magnetic field 
and it is possible for the current to depend on the sign 
of Sz. However, in cases like those considered above, 
in which this effective vector potential acts only across 
a single tunnel junction and does not produce any mag- 
netic field; the current does not depend on the sign of Sz^ 
only on its magnitude. Further, in many cases with non- 
vanishing effective field, the current will still depend only 
on the sign of Sz, as must be determined on a case- by-case 
basis; for many systems, the magnetoresistance curve is 
symmetric in the magnetization, for example [17]. 

Example with Modulation Let us finally give an ex- 
ample that does leads to a modulation at the Larmor fre- 
quency. Suppose the Hamiltonian is Ti.L+'Hn.+'HT+BSx, 
with Ht = toTo + t^JaiSpv^) + f^^MS'^^w^) and Ho = 
g°^BfSf + g^B^S'l, where there are two impurity spins, 
0,6, with different fl'-factors, coupled to the current. 
Then, choosing w to lie along the ^;-direction, the cur- 
rent modulation is 

{mm = liuifmfm^), m 

where we define n = t^S"^ + t^S''^. Then, if the two 
spins have different precession frequencies oj]^'^, the two 
spins will "beat" together and peaks will be observed 
at ±2(jj£, ±2(jj^, ±uj1 ± oj'l. If one spin precesses slowly, 
« 0, then this gives a peak at wf;. Alternately, the 
same result can be obtained by tunneling spin-polarized 
electrons, again giving rise to a modulation at 0;^. 

Conclusion — In conclusion, we have shown that the 
assumption of time reversal symmetry severely constrains 
the ability to measure a spin via a tunneling current. 
For weak magnetic fields, where we can neglect possible 
Aharonov-Bohm phases, the modulation is observed at 
twice the Larmor frequency, rather than at the Larmor 
frequency. This result is fairly general for systems with 
a single tunnel junction; for systems with a quantum dot 
in the sequential tunneling regime, or other systems in 
which the tunneling current is largely independent of the 
voltage, a peak in the current-current correlator at the 
Larmor frequency is possible. While we have derived this 
result for a specific Hamiltonian (2), it holds for a more 
general Hamiltonian, if we take H = Hl + 'Hr -\- Ht + 
BSx, with TYt given by Eq. (1). 

Given that such a modulation is observed at the Lar- 
mor frequency, and assuming the Aharonov-Bohm phases 
are negligible, then either the system is in a sequential 
tunneling regime or else there are multiple spins present, 
enabling a beating of the precession frequencies of the 



two spins. However, in the latter case there should be 
a current modulation at twice the Larmor frequency, at 
least for impurity spins greater than 1/2, and thus an 
experimental measurement of that peak is suggested. 
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FIG. 1. (a) Example diagram showing modulation at 2ojl. 
(b) Diagram which vanishes. 
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